function M2U1 = M2U1(c,gamma,cf)

M2U1=-gamma*cf*c^(-gamma-1);


end